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ABSTRACT 

A method for performing probability calculations for fuel clad surface temperatures 
of nuclear reactors is presented. The method properly treats factors which vary sys- 
tematically over the whole core, factors which vary over a plate, and those that vary at 
a spot. One method of choosing spot size for plate-type reactors is shown and a more 
general method is discussed. The method is not restricted to normal distributions. The 
calculated value for probability of success, that is, surface temperature not exceeding 
coolant saturation temperature, is identified as strictly applicable only to the instant of 
reactor startup. Also the effect of continued operation is discussed. It is observed that 
the probability of success may be more nearly equal to values obtained by the Determin- 
istic Method than to those obtained by the Statistical Method. 
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A NEW METHOD FOR TREATMENT OF UNCERTAINTIES IN NUCLEAR 
REACTOR HEAT TRANSFER CALCULATIONS 
by Harry J. Reilly, John D. Hansel I, * and George L Heath* 

Lewis Research Center 

SUMMARY 

A method of performing probability calculations for nuclear reactor surface tem- 
peratures has been devised. The method gives consideration to the fact that some uncer- 
tainties may vary systematically over all the reactor or over some parts of it. The 
method does not depend on assumption of any particular form of the probability distri- 
butions . 

The method was used to do an example calculation for an MTR-type test reactor with 
plate-type fuel elements. It was shown that the calculated probability of failure, that is, 
that surface temperature exceeds coolant saturation temperature, lies closer to values 
obtained from the Deterministic Method than to values from the Statistical Method. 

The calculated probability value was identified as the probability of success at the 
instant of reactor startup. It was observed that the probability of success for continued 
operation might not be the same as the value for startup. 

The method gives an improved representation of the probability problem for reactor 
surface temperatures. However, there is still much to be learned about the various 
important distribution functions. In the present situation of inadequate knowledge of 
behavior and distributions of uncertainty factors, all such probability calculations must 
be regarded as providing only a rough approximation to the true probability of success 
for a reactor. 


INTRODUCTION 

The designer of a high power nuclear reactor has to consider not only all the calcu- 
lated most probable temperatures and heat fluxes, but also the effects of uncertainties in 
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all the important variables. That is, the designer must allow for the fact that there are 
errors in heat transfer correlations, fuel plate loadings, reactor power level measure- 
ments, and the like. This results in the reactor being designed to operate with nominal 
or most probable values considerably less than the safe limits for those values. Of 
course, this constitutes a penalty, since it results in a reactor which operates at con- 
siderably lower power and flux levels than could be allowed if no uncertainties existed. 

During the last decade, a number of articles and reports have appeared on the 
treatment of such uncertainties by statistical methods. Nearly all of them treat the 
problem of determining the probability that the reactor fuel temperature, either clad 
surface or fuel centerline, will not exceed a limiting value, such as the saturation tem- 
perature of the cooling water. Most of these articles show the application of conventional 
statistical equations to the calculation of probability of safety for the hot spot in the core. 
These articles are reviewed in detail by Fenech and Gueron (ref. 1). 

Fenech and Gueron developed a Synthesis Method of calculation (ref. 1) . It treats 
the case of a reactor which has tubular cladding with fuel in the form of pellets. This 
type of fuel allows one to treat each pellet and its associated length of clad as a "char- 
acteristic length", where each characteristic length is independent of the others. Then 
the probability of failure for each pellet in the core is calculated. All the probabilities 
are multiplied together to get the probability of failure for the whole core, that is, 

n 

P (failure of core) = 1 - JJ [1 - P (failure of pellet n)] for n pellets (1) 

i=l 

The method is cited as having the following advantages: 

(1) It gives a realistic choice of the smallest region, or spot, that is independent of 
the other regions in the core. 

(2) It takes into account the whole core rather than just the nominally hottest spot. 

(3) The calculations, though laborious, can be done by hand, and are therefore 
useful even in the early design stages. 

(4) It does not require assumption of any particular probability distribution for the 
maximum surface temperature. 

The choice of a pellet for a characteristic length may be correct for a tubular fuel 
form. But this simple a choice is not possible for a plate-type fuel form. Also, the 
synthesis method fails to recognize that not all uncertainties are purely local in nature. 
Some quantities, such as reactor power, have error distributions which affect the 
whole core. Others, such as fuel plate loading, vary over a whole fuel plate. A valid 
treatment would have to take this into account. 

This article will show a method for performing probability calculations for fuel 
clad surface temperatures of a nuclear reactor. The method will properly treat varia- 
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bles which vary over the whole core, over one plate, and at one spot. A method of 
choosing spot size for plate -type elements will be shown, and a more general method 
discussed. All the methods are chosen for analysis using a small amount of computer 
time rather than relying on hand calculations . 

Fenech and Gueron identified existing methods as the Statistical Method, Determin- 
istic Method, Spot Method, and Synthesis Method. In this report, we preserve this 
nomenclature and refer to our proposed method as the "Probabilistic Method” for the 
sake of clarity. 

In what follows, the method of combination of probabilities to obtain the probability 
of success for the whole core is shown. Then the "spot" size derivation is given, 
followed by its use to obtain probability distributions for the "cell. " The uncertainty 
factors which can exist for an MTR-type reactor are given and an example calculation 
using the proposed method is shown. This method is compared with the Deterministic 
Method and with the Statistical Method. Finally, some observations are made about the 
significance of the values that are obtained by probability calculations of this nature. 


h 

I 0 (mr) 

^(mr) 

K Q (mr) 

Kj(mr) 

k 

m 

P(A) 

P(A|B) 

R 

RP 


SYMBOLS 

plate thickness 
constants of integration 
radial distribution function 
heat transfer coefficient 

modified Bessel function of first kind, order zero 
modified Bessel function of first kind, order one 
modified Bessel function of second kind, order zero 
modified Bessel function of second kind, order one 
thermal conductivity 
m^ = 2h/ka. 

probability of occurrence of event A 
probability of A given that B has occurred 
particular value of r 
reactor power 
radius 
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temperature 

coolant bulk temperature 
film temperature drop, T w - 
nominal value of AT^ 
maximum value of AT^ 

arbitrarily chosen member of population of AT^ 
coolant saturation temperature 
surface or "wall” temperature 
maximum surface temperature 
AT*/AT f 

fraction of heat flux at point due to heat generated within radius R of that point 
distribution function 

accumulated product of n factors 

THE COMBINATION EQUATIONS 

To better understand the method to be used here, consider a reactor composed of n 
cells, each of which operates at the same nominal temperature and each of which has the 
same distribution of uncertainty in temperature. Each cell is independent of the others. 
Now, if there is just one factor which varies in a random manner about the mean, and it 
varies independently in each cell, then 

P(success for reactor) = [P(success for cell)] n (2) 

by the multiplication rule for probability. * 

However, if this factor is one which varies in the same way over the whole core at 
once, then the cells are no longer independent, and 

P(success for reactor) = P(success for cell) (3) 

because all cells will fail at once; in fact, all will run at equal temperature. Evidently 

* "Success" means that the maximum surface temperature is less than the coolant 
saturation temperature, although this is a conservative limitation for most water-cooled 
nuclear reactors. 
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there is a difference in the effect of an uncertainty in a factor which varies over the 
whole core at once, such as reactor total power, and one which varies locally, such as 
fule plate loading. 

In the previous example, if both factors are present - that is, uncertainty in both 
reactor power and plate loadings - then the correct combination of probabilities for l 
discrete values of the reactor power is 

l n 

P(success for reactor) = ^ IpCRPj) • [P(success for cell|RPj)] 1 (4) 

i=l ^ - 

by the decomposition rule for probabilities. 

These relations are now extended to the real reactor. Suppose that the probability 
distribution for the "cell” is known, where the cell is some portion of the length of a 
fuel plate and each cell is independent of the others . Then 

n 

P(success for a fuel plate) = ]T P(success for cell i) (5) 

i=l 

There are many fuel plates in a reactor and they have a variety of neutron flux 
levels and fuel loadings. Thus the probabilities of success of fuel plates differ even at 
the same reactor power. Then 


m 

P(success for all fuel plates | RP) = J]~ P(success for fuel plate j | RP) 

j=l 


for m plates 


( 6 ) 


In some cases this can be considerably simplified. For example, if the vertical power 
distributions in all the plates are similar in shape, differing only in magnitude, then the 
radial distribution function, including all the effects that vary between plates, can be 
described by a probability distribution F r> Let each plate share the same probability 
distribution of F r , so that all plates have equal probability of success. Then 

P(success of all plates | RP) = 

for m plates per reactor, p possible values of F f (7) 


12 P(F„ ) • P(success of one plate |'F_ ) 
j=l r J J 
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Note that this involves no assumptions about the form of the various distributions other 
than that they are independent over either cell, fuel plate, or core. 


SIZE OF A "SPOT" 

The purpose of the following derivation is to provide a more logical choice of spot 
size than used in previous work. 

Figure 1 shows a cylindrical section cut through a fuel plate, with the axis of the 
cylinder normal to the plane of the plate. The plate is considered to be homogenized 



Figure 1. - Fuel plate in vicinity of "spot". 


over its thickness. A "spot" is bounded by that radius about a point such that some 
large fraction (say a) of the heat flux at the point is due to the heat source within the 
radius. The differential equation describing the temperature distribution radially from 
the axis of the spot is 

V r, y 1 

where 


T = T (plate) - T (coolant) 


The general solution is 


T = CjK 0 (mr) + C 2 lQ(mr) + 


q MT a 


2h 


( 10 ) 
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where 


2 _ 2h 
m = — 

ka 


Now consider the two -region problem where the heat generation is finite from 0 to R but 
is zero from R to ®°. Using the following boundary conditions: 


"N 


arj 

3r 


= 0 at r = 0 


3T 2 
3r 3r 


Tj = T 2 at r = R j 


( 11 ) 


gives 

Mt 

Tl = C 2 I 0 (mr) + a-i 


( 12 ) 


where 


Mt 

^-iK.(mR) 
2h 1 


K^raRjI^mR) + K^mRjI^mR) 


If the spot radius R were extended to infinity, the terms involving the partial differen- 
tials in equation (9) would vanish, giving 


T 


t f t 

= q a 

2h 


(13) 


It follows that 


K,(raR) 

a = — — =1 

T 1 (0) I <» K 0 (mR)I 1 (mR) +.Kj(mlQl 0 (mlQ 


( 14 ) 
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where T ^(0) | is the temperature at the axis when the spot radius is R and T ^ (0) | ^ 
is the temperature at the axis when the spot radius is infinite . For a typical MTR-type 
fuel element, 


k = 1.75 W/(cm)(K) 
a = 0. 15 cm 
h = 5. 0 W/(cm 2 )(K) 
m = 6. 0 cm - * 

In figure 2, a is given as a function of R and of mR. For example, a equals 0.9 when 
the radius is 0. 53 centimeter. The diameter of the spot is about 1 centimeter. If a 
60-centimeter-long by 6-centimeter-wide plate is divided into 12 lengths, there are about 
30 spots per length. 



mR, dimensionless 


Figure 2. - Fraction of heat flux at a point due to heat 
generated within radius R of that point as a function 
of R. (MTR-type fuel elements.) 
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THE "CELL" 


Previously the "cell” was mentioned as being some portion of the length of a plate. 

In fact, we arbitrarily define it to be some portion of the length of a plate. It extends the 

width of the plate and is long enough so that there are a conveniently small number of 

cells per plate. It also is small enough that the nominal heat flux and water temperature 

2 

are nearly constant in any one cell. Thus the cell is a convenient concept to use in the 
calculations, and since it is made up of a known number of "spots” which have known 
temperature probability distributions, the distribution of maximum temperature of the 
cell can be calculated. 

For example, if there are n spots in the cell, each independent of the others, and 
each having the same probability distribution of film temperature drop AT|, then 

P(AT f in a cell < ATj) = [p(AT f at a spot < ATj )J n (15) 

Suppose the distributions for the spots are given as 4»(t) . Then 


P(AT f at a spot < AT*) = J* $(t) 


dt 


(15) 


where x = AT^/AT f and <& is the distribution of AT f /AT f for the spot. Then 


P(AT f in a cell < AT|) = 


■I / 1 


$(t) dt 


n 


(17) 


gives a probability integral table for the cell. 

There is still an arbitrariness in the result, arising from the arbitrary choice of 
a which defines the spot size . Here Monte Carlo methods appear to have a great 
advantage. That is, construct a two-dimensional finite difference heat transfer model 
for the cell. Then prepare a number of problems by random sampling from the prob- 
ability distributions within the cell. In this way, the probability distribution of the 
maximum film temperature drop factor in the cell can be calculated. Then use the 
combination equations shown here to build up the probabilities of success for the fuel 
plates and the reactor core. This work did not include a Monte Carlo calculation, be- 

2 

There may actually be some local flux peaking in the cell. But since the spots are 
independent, not all spots in the cell have to come from the same physical location. The 
cell just accounts for a certain number of spots, which are usually close together. 
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cause the statistical data are not good enough to warrant it; the spot size derivation 
shown was adequate to do the comparisons and draw the relevant conclusions. 


EXAMPLE CALCULATIONS AND COMPARISONS 

An example calculation was done using the operating parameters and the uncertain- 
ties of a typical light-water test reactor at startup. The operating parameters and the 
uncertainty factors and their values are listed in table I. 


TABLE I. - NOMINAL OPERATING PARAMETERS AND 
UNCERTAINTY FACTORS 


(a) Nominal operating parameters at startup 


Power, MW(th) 

Inlet temperature, K 

Coolant saturation temperature (at core outlet), K 


60 

330 

440 


(b) Uncertainty factors 


Region 

of 

inde- 

pendence 

Factor 

Fractional 

standard 

deviation 

percent 

Distribution a 

Spot 

Local fuel density x thickness 

5 

Normal 


Local heat transfer coefficient 

10 

Normal 


Flux peak at bottom end of fuel 

4.5 

Normal 

Plate 

Coolant velocity 

6 

Normal 


Power distribution 


Rectangular, 
1. 0±0. 5 


Fuel per plate 


Normal 


Area per plate 


Normal 

Core 

Reactor power 

5 

Normal 


^he distributions were all arbitrarily truncated at ±3a. 


The rectangular power distribution is illustrated in figure 3. It was not assumed 
that the power distribution in the core was uniform; rather, it was assumed that a given 
plate had equal probability of having any power between 0. 5 and 1.5 times the average 
power. Figure 2 also shows the normal distribution of the product of the coolant velocity, 
fuel per plate, and area per plate factors. (That distribution was also arbitrarily trun- 
cated at ±3o.) Finally, figure 3 shows the distribution of the product of the factor having 
normal distribution and the factor having rectangular distribution. The result is the 
F r distribution referred to earlier . 

Figure 4 shows the probability integral of AT f /AT f for the spot and for the cell 
when there are 30 spots in the cell. The two curves are related by equation (17). Notice 
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Value of factor divided by its nominal value, A 

Figure 3. - Probability distributions for factors which vary in- 
dependently over fuel plates. 



1.0 1.06' 1.12 1.18 1.24 1.30 1.36 

N 

Figure 4. - Probability that ATf/ATj is less than N in spot 
and in cell (30 spots per cell!. 
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that there is negligible probability that AT f /AT f will be less than 1. 12 for the cell. No 
value can be greater than 1. 36, because the distributions are truncated at ±3o (three 
times the standard deviation) . 


Figure 5 shows the calculated nominal values of T gat , T^, and T w 


when F r = 


1. 0 or 1.8. The numbering of the cells is shown along the abscissa. The value of the 
ratio (T gat - Tk)/(T W - T^) = 5.75 at cell 4 when F r = 1.0. The probability that 
AT f /AT f exceeded 5.75 was zero from figure 4. Thus the probability that T w < T gat 
was 1.0 for this cell. The probability that T w was less than T oot was calculated for 


sat 


Radial 



Figure 5. - Temperature as a function of position in channel 
for two values of F f . (Reactor at 60 MW. I 


each cell and the values combined as in equation (5) to obtain the probability that 
T w < T gat for the whole plate. This was done for each value of F r at each reactor 
power. Figure 6 shows P(T w < T gat ) as a function of % F r at 60 megawatts. This 
result for each reactor power was used as in equation (7) to obtain figure 7 . Figure 7 
shows P(T W < T ga j.) for the whole core as a function of true reactor power. The values 
given in figure 7 were then used in equation (8) to obtain P(T W < T gat ) for the reactor, 
including the effects of uncertainty in reactor power level. 

The uncertainty factors of table I were also used in calculations by the Deterministic 
Method and by the Statistical Method. (In the Deterministic Method, the maximum possi- 
ble values are all assumed to occur at once; the +3o values are the maximum possible 
in this example. In the Statistical Method, the uncertainty in the temperature is calcu- 
lated only for the nominally hottest spot.) 

The results of calculations by the Deterministic Method and by the Statistical Method 
are compared with the result of the new Probabilistic Method given here in table n. The 
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1.0 



0L- 

1.0 


1.2 


_J I 

1.4 1.6 

Radial factor, F r 


Figure 6. - Probability that T w <T sa ( for one 
fuel plate as a function of radial factor F r . 
(Reactor at 60 MW. ) 



Figure 7. - P (success) for core as a function of true reactor 
power. 
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TABLE n. - COMPARISON OF RESULTS USING 


DIFFERENT METHODS 


Method 

Results 

P(success) 

Deterministic 

Statistical 

Probabilistic 

h varying locally 
h varying over whole 
core systematically 

T w = 485 K, T sat = 440 K 
Most probable T w = 409 K 

0.9995 

.75 

.92 




Probabilistic Method gives a probability of success of 0. 75. The Deterministic Method 
gives a value of 485 K, which is greater than T ga ^.. According to the Deterministic 
Method, there is zero probability that T w exceeds 485 K and unknown probabilities for 
all lesser values. 

The Statistical Method gives a most probable T w of 409 K with a standard deviation 
of 9.7 K. This has a normal distribution such that the probability of success is 0.9995. 
Thus P(failure) = 1 - 0.9995 = 0.0005 compared to 1 - 0.75 = 0.25 for the Probabilistic 
Method. The results differ by a factor of 500. The Deterministic Method in effect gives 
a P (failure) of 1.0. This result for P (failure) is greater than that of the Probabilistic 
Method by a factor of 4.0. Thus the Deterministic Method, which is simpler and more 
conservative than the Statistical Method, also is in better agreement with the Probabil- 
istic Method, at least for this example . 

To illustrate the effect of core-wise variation instead of local variation, the Prob- 
abilistic Method calculations were redone with one change. The heat transfer coefficient 
h was assumed to vary systematically over the whole core instead of locally. The 
result was a P(success) =0.92 for the core. The probability of failure is thus about 
three times smaller if h varies core wise rather than locally. This demonstrates the 
importance of knowledge of the manner of variation of each factor and the importance of 
proper combination of the different uncertainty factors. 


EFFECT OF CONTINUED OPERATION 

The probability distribution of maximum surface temperature is a function of time. 
In an MTR-type reactor, for each core loading, the distributions of fuel and the other 
variables in the analysis stay about the same as they are at the start of the operating 
cycle. The neutron flux peaking decreases during the cycle because of control rod 
withdrawal, so the maximum surface temperature T w tends to decrease with time. 
This is illustrated in figure 8. 
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The value 0.75 that we obtained for P(success) of the core is equal to 



( 18 ) 


in figure 8. There is 0.75 probability of starting the cycle in the ’’safe" zone, say at 

A 

point P. If conditions stay the same, T should progress along path A. If some things 

A ” 

change in a random manner, T w will follow a ’’random walk” from left to right in fig- 
ure 8, with a randomly chosen step toward larger or smaller T w during each step in 
time. But an accurate construction of this random walk would require knowledge of how 
rapidly and extensively the variables could change with time, and such knowledge does 
not now exist. 

For the example, it seems reasonable to assume that the calculated value of 
P(success) equals that for one cycle of operation. But each cycle may be a repeat of the 
trial because a new core loading is used for every cycle. If so, P(success) for 10 cycles 
is (0.75)^ or about 0.06. 

There are two conclusions to draw. First, the value for P(success) calculated by 
this and most other existing probabilistic or statistical methods is identified as appro- 
priate to the instant of startup. Secondly, the P(success) for any continuing operation 
may be less than that for the instant of startup. This aspect of the problem must be 
considered by the designer when using such methods. 



Figure 8. - Probability distribution of maximum surface temperature T w 
as a function of time t. 
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CONCLUSIONS 


A method of performing probability calculations for nuclear reactor surface tem- 
peratures has been devised. The method gives consideration to the fact that some un- 
certainties may vary systematically over all the reactor or over some parts of it. The 
method does not depend on assumption of any particular form of the probability distri- 
butions. 

The method was used to do an example calculation for an MTR-type test reactor 
with plate-type fuel elements. It was shown that the calculated probability of failure, 
that is, that surface temperature exceeds coolant saturation temperature, lies closer to 
values obtained from the Deterministic Method than to values from the Statistical Method. 

The calculated probability value was identified as the probability of success at the 
instant of reactor startup. It was observed that the probability of success for continued 
operation might not be the same as the value for startup. 

The method gives an improved representation of the probability problem for reactor 
surface temperatures. However, there is still much to be learned about the various 
important distribution functions . In the present situation of inadequate knowledge of 
behavior and distributions of uncertainty factors, all such probability calculations must 
be regarded as providing only a rough approximation to the true probability of success 
for a reactor. 
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